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Abstract 

SQUID magnetization measurements in oriented pow- 
ders of Yi_a;Ca2;Ba2Cu30y, with x ranging from to 0.2, for y « 6.1 and 
y ~ 6.97, have been performed in order to study the doping dependence of 
the fluctuating diamagnetism above the superconducting transition tempera- 
ture Tc- While for optimally doped compounds the diamagnetic susceptibility 
and the magnetization curves —Mfi{T = const) vs. H are rather weU jus- 
tified on the basis of an anisotropic Ginzburg-Landau (GL) functional, in 
underdoped and overdoped regimes an anomalous diamagnetism is observed, 
with a large enhancement with respect to the GL scenario. Furthermore the 
shape of magnetization curves differs strongly from the one derived in that 
scheme. The anomalies are discussed in terms of phase fluctuations of the 
order parameter in a layered system of vortices and in the assumption of 
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charge inhomogeneities inducing local, non percolating, superconducting re- 
gions with Tc°'^^ higher than the resistive transition temperature Tc. The 
susceptibihty displays activated temperature behavior, a mark characteristic 
of the vortex-antivortex description, while history dependent magnetization, 
with relaxation after zero-field cooling, is consistent with the hypothesis of 
superconducting droplets in the normal state. Thus the theoretical picture 
consistently accounts for most experimental findings. 



PACS: 74.40.+k, 74.70.Vy,74.25.Ha 
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I. INTRODUCTION 



A variety of experiments^ points out that the small coherence length, reduced carrier 
density, high transition temperature Tc and marked anisotropy of cuprate superconductors 
cause strong enhancement of superconducting fluctuations (SF). In contrast to conventional 
superconductors, in cuprates the transition region is considerably smeared by SF which can 
be detected in a wide temperature range, up to 10-15 K. The formation of the fluctuation 
Cooper pairs above results in the appearance of a Langevin-type diamagnetic contribution 
to the magnetization —Mfi{T,H), existing side by side to the paramagnetic contribution 
from fermionic carriers. 

Since the size of fluctuating pairs ^ (T) grows when T approaches the transition temper- 
ature Tc, Mfi{T,H) should diverge near the transition for any small fixed magnetic field, 
being equal zero for if = 0. On the other hand it is evident that very strong magnetic fields, 
comparable to Hc2{0), must suppress SF. Therefore the isothermal magnetization curve 
Mfi (T = const, H) has to exhibit an upturn. This upturn can be quantitatively described 
in the framework of the exactly solvable, for any magnetic field, zero-dimensional modeP 
(superconducting granula with the size <S ^{T)) or by means of cumbersome microscopic 
treatment accounting for the short-wavelength fluctuation contribution in the ?>D case.^ 

The experiments on conventional BCS superconductors show that the magnetization is 
quenched for fields as low as ~ 10~^i?c2(0) (see Ref.^). The value of the upturn field in 
the magnetization curves can be considered inversely proportional to the coherence length. ^'^ 
This explains why in optimally doped high-temperature superconductors, the Ginzburg- 
Landau (GL) picture works pretty well. Here the coherence length is so short that the 
quenching of fluctuating magnetization on increasing the magnetic fleld has not yet been 
observed. 

The fluctuating magnetization of layered superconductors in the vicinity of the transition 
temperature and for iic2(0), when the contribution of short-wavelength fluctuations'^ is 
negligible, can be theoretically described^"'^ in the framework of the GL scheme with the 
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Lawrence-Doniach Hamiltonian.^"^'^ The fluctuating diamagnetism (FD) turns out to be a 
complicated nonlinear function of temperature and magnetic field and cannot be factorized 
on these variables. An important role in FD is played by the degree of anisotropy of the 
electronic spectrum. All these aspect of FD have been found to occur in optimally doped 
YBCO.^~^^ Also scaling arguments^^ were found^^'^^ rather well obeyed in this compound. 

In underdoped YBCO, instead, marked deviations from the behaviour expected in the 
framework of GL approaches have been detected. A first qualitative claim in this regard 
goes back to Kanoda et al./^ who noticed that in oxygen deficient YBCO the FD in small 
fields was enhanced. Later on, novel features of FD in underdoped compounds have been 
reported ^'^6-20 particular, in underdoped YBCO at Tc ~ 63^ marked enhancement 
of the susceptibihty for fixed (T — in a field of 0.02 Tesla was detected above Tc(0), 
with magnetization curves strongly different from the ones in optimally doped YBCO.^^'^^ 
Magnetization curves have been subsequently reported in underdoped Lai,gSrQ,iCu04 
(LASCO). For the moment we only mention that the magnetization curves reported by 
Carballeira et al.^^ in underdoped LASCO, although indicating field-affected fluctuation- 
induced diamagnetism, do not exhibit the upturn with the magnetic fleld as the ones in 
underdoped YBCO that we will discuss later on. Finally, recent magnetization data^° as a 
function of temperature in YBa2CusOQ,5 single crystal (with transition temperature in zero 
fleld T'c(O) = indicate SF obeying to 2D scahng conditions for HZl Tesla and turning 

to 3D scaling for smaller flelds. 

Qualitative justiflcations of the anomalous diamagnetism in underdoped YBCO have 
been tried^^'^^, essentially based on the idea of charge inhomogeneities leading to non- 
percolating superconducting "drops" or on the extension of the theory by Ovchinnikov et 
al.^^, where the anomalous diamagnetism is related to regions having local T^'s higher than 
the resistive transition temperature. The flrst theoretical study speciflcally aimed at the 
description of FD in underdoped YBCO was undertaken by Sewer and Beck In the 
framework of the Lawrence-Doniach model, these authors justify the temperature and field 
dependences of the magnetic susceptibility by taking into account the phase fiuctuations of 
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the order parameter, thus arriving to a layered XY-model for a hquid of vortices. 

In this paper we address the problem of the fluctuating diamagnetism in the 
Yi_j;Caa;Ba2Cu30y family and of its dependence from the number of holes, by reporting 
SQUID magnetization measurements in a series of samples. Preliminary results on over- 
doped compounds have been presented to a meeting and published elsewhere. Since some 
differences in the magnetic behaviour of chain-ordered and chain-disordered YBCO have 
been noticed^^, we also attained the underdoped regime by means of Cc?^ for sub- 
stitution in ideally chain-empty YBa2Cu30Q, while for the overdoped regime the same 
heterovalent substitution was performed in chain-full YBa2Cu307. 

The paper is organized as follows. In Sect. II experimental details and the majority of 
the experimental results are reported. The analysis of the data (Sect. Ill) is first tentatively 
carried out on the basis of an anisotropic free energy GL functional, within the Gaussian 
approximation. The inapplicability of such an approach for non-optimally doped compounds 
is stressed. Then the theory for phase fluctuations of the order parameter in layered liq- 
uid of vortices is revised, to properly take into account terms neglected in the previous 
formulation.^^ In particular non-reversibility and relaxation effects of the magnetization are 
argued to support the picture of non-percolating, locally superconducting droplets above the 
resistive transition temperature, that we interpret as phase fluctuations of a non-zero order 
parameter below the local irreversibility temperature. Thus a comprehensive description of 
FD in the Yi_j;Caj;Ba2Cu30j/ family is obtained, as it is summarized in Sect. IV. 

II. EXPERIMENTAL DETAILS AND RESULTS 

The samples of chemical composition Yi_xCaxBa2CusOy were prepared by solid state 
reactions of oxides and carbonates in flowing oxygen at 1000 K for about 100 hours. X-Ray 
diffractometry was used to check the presence of a single phase. The oxygen stoichiometry 
was first estimated by thermogravimetry and energy dispersive spectrometry. The samples 
were then oxygenated close to y = 7 by annealing in oxygen atmosphere (25 Atm) at 450 K 
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for about 100 hours or deoxygenated as much as possible close to y = 6, for about 100 hours 
in vacuum. The final oxygen content turned out y — 6.97 ± 0.02 for overdoped YBCO and 
7/ = 6.10 ± 0.05 for the underdoped samples, estimated with loss of mass measures. Before 
the measurements, the samples stay at room temperature for about 1 week. The resistive 
transition in samples of the same batch appeared very sharp, with moderate evidence of 
paraconductivity in a temperature range of 5-10 K above the transition. After mixing the 
samples with epoxy resin, they were oriented by hardening in a strong magnetic field ( 9 
Tesla). The orientation was tested by comparing the diamagnetic susceptibility for H//c 
with the one for H in the ab plane, where practically no enhancement of M/; was noted to 
occur. For samples used in previous works the orientation was also tested by means of 
the 6^Cu NMR line (see ref. 18). 

Magnetization measurements have been carried out in the oriented powders by means 
of a Quantum Design MPMS-XL7 SQUID magnetometer. Measurements were performed 
also in optimally doped YBCO in order to prove that the results in oriented powders do 
not significantly differ from the ones in single crystals. The data already obtained by other 
authors^'^'^^'^^"^^ were confirmed. In the following Section we will recall a few results of the 
studies in optimally doped compounds, when it is required for the comparison of our data 
in strongly underdoped or overdoped YBCO:Ca. 

The transition temperatures TjyE. = 0) = T'c(O) have been estimated from the mag- 
netization curves vs. T at small fields (20 Oersted), by extrapolating at M=0 the linear 
behavior of x occurring below Tc, as shown in the insets in Fig.l. The values of Tc(0) 
are collected in Table I, where the numbers of holes n/j, as evaluated from the expression 
(r^) — 1 ~ 82.6(n/i — 0.16)^, giving the parabolic behaviour^^ of the phase diagram T — x, 
are also reported. It is noted that because of the enhanced fluctuating diamagnetism some 
uncertainty in the estimate of Tc(0) is present, particularly in strongly underdoped sam- 
ples. This uncertainty does not affect the discussion given later on about the anomalous 
diamagnetism, which is detected in a temperature region well above Tc(0). 

Magnetization measurements at constant field H have been performed as a function 
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of temperature, with H || c. In general two contributions to the magnetization M were 
observed: a Pauh-hke, positive term Mp, almost T-independent or only slightly increasing 
with decreasing temperature in the range AT from 200K down to about 100 K and a negative 
diamagnetic contribution —Mfi arising on approaching Tc . This latter contribution was 
extracted by subtracting from M the value obtained by extrapolating for T — > the 
curve Mp vs. T in AT, where Mfi is practically zero. Thus the possible slight temperature 
dependence of Mp around Tc was neglected in comparison to the much stronger diamagnetic 
term. 

Typical magnetization curves M{H = const, T)for overdoped and underdoped samples 
are reported in Fig.2. The enhancement of FD, in both regimes, is evidenced. 

In Fig. 3 some isothermal magnetization curves —Mfi{T — const) vs H, obtained by 
cooling in zero magnetic field (ZFC) down to a certain temperature above Tc(0), are shown. 
In Fig 3c also a few data for field cooled (FC) magnetization, to be discussed later on, are 
reported. 

III. ANALYSIS OF THE DATA, FURTHER RESULTS AND THE 

THEORETICAL PICTURE 

A. GL anisotropic free energy functional 

The generalization of the GL functional for layered superconductors [Lawrence-Doniach 
(LD) functional ^'^^] in a perpendicular magnetic field can be written 

(-"-^-")*'r 

(1) 

where \E'; is the order parameter of the /— th superconducting layer and the phenomenological 

constant is proportional to the Josephson coupling between adjacent planes and a — 
T — Tc 

ao{ — — — -) = ckqE. The gauge — is chosen in Eq.(l). In the vicinity of Tc the LD 

c 

7 



LD 



167r2 



m 



functional is reduced to the GL one with the effective mass M — (4J's^)~"^ along c-direction, 
where s is the inter-layer spacing. In the GL region the fourth order term in (1) is omitted 
and the standard procedure^'"^ to derive the fluctuation part of the free energy yields 

F(e,iJ)-F(e,0)= (2) 
TVkB , r , f^^^ , , (2n + 1 + 2x)h + r/2(l - cos z) + e 



27rsC, 



, r J V- f ' , , 2n + l + 2x /i + r/2 1-cos^ -+ 

„ h dzy / dxm — -— — r 

I, (2n+l)/i + r/2(l-cosz)+e 



where the c-axis is along the z direction, r — and h — jf^^^j- 

By means of numerical derivation of Eq.2 with respect to the field one obtains the 
fluctuating magnetization Mfi vs. H. As shown in the inset of Fig.4, the magnetization 
curves in optimally doped YBCO are satisfactorily fitted by Mfi derived in this way and 
evidence how the 3D scenario of SF is obeyed on approaching Tc, with a crossover from linear 
to non-linear field dependence occurring a few degrees above the transition. Correspondingly, 
the scaling arguments for 3D anisotropic systems hold and Mfi/H^/'^ vs. T cross at Tc(0) ~ 
927^, as already observed^^. 

In contrast to optimally doped YBCO, the magnetization curves for underdoped and 
overdoped compounds depart in a dramatic way from the ones expected on the basis of 
Eq.2. In particular (see Fig.2a-3b-3c), even relatively far from Tc, while for small fields 
(H(i71inear in H, upon increasing the field the magnetization shows an upturn and then 
I — Mfi I decreases. Let us remind that in the GL weak fluctuation regime the saturation of 
the magnetization at high field has to be expected^°'^^, the superconducting coherence being 
broken for fields larger than ^/eHc2■ An estimate of the order of magnitude of the upturn field 
Hup can be done from the analysis of the 'OD' case^'^, namely for superconducting granules 
of radius d smaller than the coherence length ^(T). In this case the order parameter is 
spatially homogeneous and the exact solution of the GL model can be found and yields 



Mfi = . (3) 

It can be noticed that the most sizeable contribution to the magnetization comes from the 
fluctuations-induced SC droplets of radius of the order of ^ (T) , which imply most efficient 



screening. By assuming the condition of zero dimension for these droplets, from Eq.3 with 
d = ^{T) one derives an upturn field given by H^p ^ For { ~ lOA and e in the range 

10"^ — 10~^ , Hup is expected to be in the range of ten Tesla. 

Thus the magnetization curves in Fig.s 3 can hardly be ascribed to the breakdown of 
the GL approach of the type commonly observed in BCS superconductors^. In other words, 
a description of FD based on the GL functional in principle should be suitable in YBCO 
compounds for fields smaller than several Tesla, particularly not too close to Tc, as in fact 
it is observed in optimally doped YBCO.^'^'^^'^^'^''' 

B. Phase fluctuations and superconducting droplets above Tc : a theoretical picture 

Then one has to look for other explanations. As already mentioned, a recent theory 
has been developed by Sewer and Beck^^ at the aim to justify the unusual magnetization 
curves detected in underdoped YBCO^^. The theory assumes a frozen amplitude of the order 
parameter while phase fluctuations are taken into account. As a consequence one has to 
deal both with thermally activated vortex loops and field induced vortex lines. Two major 
conclusions of general character can be outlined. For small field the temperature dependence 
of the susceptibility is controlled by the vortex loops density n^, for which 

= noexp[-Eo/kT] (4) 

according to the XY model. For strong fields, instead, the vortex hne elements dominate, the 
vortex correlations between different layers become relevant and Mji only slightly increases 
with H and finally it flattens. No upturn field is predicted, at least for i7 << Hc2 ■ 

In Fig. 5 the data for x, defined as {—Mfi/H), are shown to obey rather well to Eq.4, 
in correspondence to Eq ~ 9A0K, in agreement with calculations yielding for the activation 
energy values around 10 Tc (see Ref.^^ and references therein). Eq turns out to depend 
only little from doping, being slightly field dependent. It is necessary to mention that the 
temperature dependence of the susceptibility above Tc(0) differs from that one measured 
below Tc, as shown in Fig. 5. 
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However the magnetization curves, as the ones reported in Figs. 3 and 4, cannot be 
accounted for by a theory which does not include an upturn with the field. Furthermore 
similar effects are found also in overdoped samples (see Figs. 2a and 3a). 

Thus we are going to consider the second aspect possibly leading to an anomalous 
diamagnetism, the one related to charge inhomogeneities causing regions where the hole 
density is different from the average. Evidences of inhomogeneous structure of cuprates 
have been found by means of neutron and electron diffraction^°~^^, as related to stripes 
and lattice effects or to local variation in the oxygen concentration, particularly near grain 
boundaries. Intrinsic inhomogeneities, with spatially dependent critical temperature have 
been considered as possible cause of pseudogap phenomena^^. A theory for high temper- 
ature superconductivity and of the pseudogap temperature dependence based on inhomo- 
geneous charge distribution with site-dependent transition temperature has been recently 
formulated.^^ In particular, Ovchinnikov et al.^^ have considered the anomalous diamag- 
netism above T^. induced by non-uniform distribution of magnetic impurities, depressing 
but leaving "islands" which become superconductors above the resistive transition temper- 
ature. An anomalous large diamagnetic moment results above Tc and in this way the strong 
diamagnetic susceptibility observed in overdoped Tl-based cuprates^^ could be explained. It 
should be stressed, however, that in this description^^ the magnetization is linear in the field, 
since the condition of small field is implicitly assumed. Direct evidence of inhomogeneous 
magnetic domains showing diamagnetic activity above Tp has been obtained by Iguchi et 
36 scanning SQUID microscopy in underdoped LASCO. Regions of few tens of yum, 
precursors of bulk superconductivity have been imaged in this remarkable work.^^ 

In the light of the experimental findings and of the theoretical supports outlined above we 
consider now as a source of diamagnetism above T^. the presence of locally superconducting 
droplets. From the volume susceptibility, let us say at Tg— SiC (see Fig.l), one deduces that a 
few percents of the total material being superconductor above the resistive transition, could 
actually justify the screening effects observed as FD. A test of this hypothesis is obtained 
from the search of magnetic-history dependent effects. It is known, in fact, that in YBCO 
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the irreversibility temperature is not far from Tc and therefore if the anomalous FD has to 
be attributed to locally SC droplets then one should detect differences between ZFC and 
field-cooled (FC) magnetization. In Fig. 6 magnetization curves after zero field cooling 
and the correspondent values of Mfi obtained at the same temperature after cooling in the 
presence of a given magnetic field, are compared. Furthermore relaxation effects have been 
observed. In Fig. 7 it is shown how the negative magnetization depends on time, displaying 
a progressive decrease from the ZFC value towards the one measured in FC condition. The 
time constant for this relaxation process is close to the one measured in the critical state^^. It 
can be remarked than in underdoped chain-disordered YBCO (Fig. 3c) no upturn is observed 
and the ZFC and FC magnetization curves almost coincide. The explanation that will be 
supported from our theoretical picture is that the magnetization curves without hysteretic 
effects refer to superconducting droplets which are above the irreversibility temperature. 

One could suspect that the occurrence of superconducting droplets results from trivial 
chemical inhomogeneities of the samples. As described in the Section on experimental details 
many experimental checks allow us to rule out this hypothesis. Furthermore, samples grown 
with different procedures and already used by other authors, have been studied. Thus it 
is believed that the inhomogeneity does not mean the presence of macroscopic parts of 
the samples at different oxygen and/or calcium content, but it is rather intrinsic, as the 
ones evidenced in the experiments recalled above. Furthermore it should be remarked that 
the temperature dependence of the susceptibility above the bulk transition temperature is 
different from the one occurring in the superconducting state (see Fig. 5). 

In the following we are going to modify the theoretical description of Sewer and Beck^^, 
still keeping their basic idea of phase fiuctuations but taking into account the presence of 
mesoscopic "islands" with non-zero average order parameter amplitude that can be below 
or above the local irreversibility temperature. 

Let us start, as in Ref.^^, from Eq.l by evidencing the order parameter phase contribution 

J'ld [0] = / I (^"^ - ^ A") + -^^[1 - cos(^/+i - ei)] I (5) 
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where J„ = and J±^ — zTrjrih are the order parameter phase couphng constants on 

ArUe 

the plane and between planes respectively. 

In this way the occurrence of superconducting droplets below the critical temperature 

is assumed, where the order parameter phase can fluctuate producing thermal excitations 

(vortex and antivortex pairs in 2D, vortex loops in anisotropic model). The potential vector 

A|| in Eq.5 describes both the magnetic field applied parallel to the c - axis and the one 

induced by thermal fluctuations. 

By following the 2D Coulomb gas theory, at each vortex is associated an effective charge 

J- 

Qv — -\/27rX and a vortex-antivortex pair has an energy Eq — q^\n{—), playing the role of 
an activation energy and thus yielding Eq.4. In order to refer to the anisotropic 3D model 
the vortex lines (or the vertical elements of the vortex loops) are correlated along the c - 
axis for a length ns and a correction to was found selfconsistently. 

By considering, as usual, the partition function Z = J DO exp{—i3J^Lc> [0]) with j3 — 

1 dMfi dF 

, the susceptibility x = , where Mfi = ——, is obtained as the sum of three 



keT' ^ ' dH ' dH 

contributions: 

where () means the thermal average. 

In the gauge A — —yH, z being the c-axis direction, the homogeneous susceptibility is 
given by 

X =lto ^, (7) 



where 

T /9^\ 2 r 7" 1 

(8) 



^(P(q)-Q(q))-l 



In Eq.(8) P{q) derives from the term {{ ^^'^ Y) of Eq.(6) and it involves the current- 
current correlation function, as in Ref.^^ : 
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^(q) = E / ^'/^ / d'p' exp [^q (r - r')] 
1,1' 



(9) 



with N the number of layers and = Trit!^, R being the average radius of the supercon- 
ducting islands. 

The x-component of the phase gradient is 

^xen{p) = dy y-Ry(mi,h) — t{rm,h), 

where t (mi, h) ~ ±1 and R(mi, /i) labels the position of each "pancake" mi on the layer 
k. 

Three terms are obtained by the evaluation of Eq.(9): Peeio) , Paa{q) and P0a{q) (the 
two terms due to the correlation between Vxd and A,, -r give the same contribution being 
PeA{q) — -Pa6>(?) — PeA{—q))- The first one involves the positional correlation function of 
the vortex line elements. In order to calculate it, Sewer and Beck^^ introduced the static 
structure factor of a disordered vortex hquid. Because of the weak interlayer coupling 
harmonic deviations of the vortex lines (or loops) along the z direction are taken into 
account. This model can be used to describe also the vortex system in the glassy phase, 
below the irreversibility line temperature and therefore the same expression for Peg{q) is 
used here. 

The evaluation of the term Paa{q) is straightforward and one has 

The further contribution PeAi^) , appearing due to the cross correlation between S/xO 
and Alia; and disregarded in Ref.^^, cannot be neglected below the vortex lattice melting 
temperature, where irreversibility effects occur. In this case one obtains 

PeA{q) + PeA{-q) = 2—^^ ^— ^ ^ exp {-dq \l-h\) 

^ *o ^ 1,1' 



t (mi, li) COS [iqlRy (mi, li)] \ . 
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The thermal average is performed in the assumption that the vortices are uniformly 
distributed in the planes and the calculations are reported in Appendix. The expansion of 
P0a{q) in powers of qL gives 



PeA — 



27r2 rHL 



-2\ 2 



+ — TT^ 



(11) 



3 V ^0 / 45 V $0 
The function Q{q) in Eq.8, related to the third term of Eq.6, has been neglected in Ref.^^. 
It can be calculated as described in Appendix, yielding 



Q = (27r: 



1 1 2.2 , 1 

^2^2 + X 4^ 12 



(12) 



It should be noted that the first term in Q, diverging for g — > 0, exactly cancels out the 
term in the expansion of P{q) which appears from the structure factor. 

By using Eqs. (10), (11), (12) and Eq.(6) of Ref.^^, from Eq.(8) one finally obtains 



K{q) 



■Li ( ^ 

s V$o 



H 

IT* 



- 1 



+ 



(13) 



kT 



^ 1 (l + ^(#^)> _ fV(l + n) (H_\\ 
ll + 2n l + 2n ^ V^V 



540 s V^oi '^i^ 



with 5 = ^—--^ and H 



is an effective " critical" field depending on the island size. 



4 kT " " ttR^ 

To avoid unphysical divergences in the calculation of the susceptibihty from Eq.(7), the 
first term in square brackets of Eq.(13) has to be zero, giving a renormalization of q^ due to 
both the anisotropy of the system and the presence of applied magnetic field: 



(14) 



In view of the field-dependent vortex charge, the pair energy (in the limit H < H* ) 
becomes E = 5-. According to Eq.(4) the thermally-excited vortex pair density 

i^ + Hw) ) 

turns out field dependent. This field dependence, formally derived in our description, is 
significantly different from the one assumed in Ref.^^. 
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Finally the diamagnetic susceptibility is obtained in the form 



^~ s^ll + 2n riy l + 2n ^ ^ \H* J 540 s \%) \H* 

(15) 

In the limit H ^ a good agreement of the susceptibility and its temperature dependence 
with the experimental findings is again achieved. The main differences between our suscepti- 
bility in Eq.(15) and the one given in Ref.^^ consists in the presence of the factor ij^Y 
of the third, positive term. This term can give an inversion in the sign of the susceptibility 
corresponding to an upturn in the magnetization curves. This phenomenon depends on the 
dimension of the islands and x = (i.e. the occurrence of the upturn) requires R > Rq 
where Rq depends on some characteristics of the material. By choosing 7 = 6, the interlayer 

distance s = 12A, n = 2, = 2.5, which are typical values for YBCO, for T = 75.5K 

kT 

one estimates Rq ^ 50A. In this case the solutions of the Equation x = is ~ 0.06 
and by considering the experimental value H^p ^ 250 G, the effective critical field turns out 
H* ~ OAT. 

The isothermal curves can be obtained from Eq.(15) by means of numerical integration. 
The shape of the magnetization curve depends on the parameters in the susceptibility and 
by using the values quoted before, with R — 370A one derives the behaviour sketched in 
Figs. 3a for an island below the irreversibility line. The same parameters with T — 66K, 

"4^ = 1.8 and it! = lOA lead to the curve shown in Fig.Sc for the magnetization of the island 

kT 

above irreversibility. 

Finally we discuss the differences observed in the magnetization curves between chain- 
ordered and chain-disordered YBCO compounds and the relevant observation by Carballeira 
et al.^^ of magnetization curves -Mfi vs H in underdoped LASCO (Sr content a;=0.1, 
Tc(0)=27.1 K) with no upturn field. A role of the chains on favouring the nucleation of 
local superconducting droplets above Tc is conceivable. In fact, in chain-ordered compounds 
the droplets appear to have the irreversibility temperature higher than the ones in chain- 
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disordered compounds, as it is evidenced by the difference in the magnetization curves (see 
Figs. 3a-b-c). We remind that the inversion in the sign of the susceptibihty is related to 
the third term in Eq.l5 and thus to the term PgA.{q)- The amount of impurities and/or 
imperfections acting as nucleation centers might also play a role. Furthermore the degree 
of under or over-doping is also involved since a marked variation of Tc with n^ is evidently 
crucial. It is noted that in the measurements by Carballeira et al.^^ in LASCO at Tc=27.1 
K the magnetization curves show only a weak tendency to saturation while in LASCO at 
Tc ~18 K (therefore more underdoped) scanning SQUID microscopy by Iguchi et aP^ does 
evidence diamagnetic effects to associate to locally superconducting droplets. These droplets 
should imply a contribution to the magnetization curves similar to the one detected by us 
in YBCO compounds. Future research work will have to explore these interesting aspects 
and the differences until now present between LASCO and YBCO. 
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IV. CONCLUSIONS 



By means of SQUID measurements in Yi_xCaxBa2CusOy family a non conventional 
fluctuating diamagnetism has been observed in overdoped and in underdoped compounds. 
Compared to optimally doped YBCO, a large enhancement of the diamagnetic susceptibil- 
ity occurs and no anisotropic GL functional or scaling arguments can justify the isothermal 
magnetization curves. The recent theory^^ for phase fluctuations of the order parameter 
in a layered liquid of vortices has been revised and it appears to justify some aspects of 
the anomalous FD in non-optimally doped YBCO, particularly the "precritical" tempera- 
ture activated behaviour of the susceptibility in the limit of zero field. Other experimental 
observations, noticeably the upturn in the field dependence of the isothermal fluctuating 
magnetization and history-dependent effects, indicate the role of mesoscopic charge inho- 
mogeneities in inducing local, non-percolating, superconducting "droplets". On the basis 
of both types of experimental findings we have extended the theory of phase fluctuations 
in the presence of non-zero order parameter. The terms leading to a novel and relevant 
dependence of the fluctuating magnetization from the magnetic field were included in the 
scheme. The field-related corrections are different when the superconducting droplets are 
below or above the local irreversibility temperature. In this way most of the experimental 
findings have been justified. 
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VII. APPENDIX 



To derive Eqs.(ll) and (12) the following thermal average must be calculated: 

^ t (m, /) cos iqlR, (m, I)) =Y, (^M) " 5^ < ^'(^) > +^(^') (^1) 

^ m,l / m m 

Indicating with A^_|_ {N_) the number of vortex line elements parallel (antiparallel) to the 
field the first term gives 

00 

The sum in the second term can be split in two parts which separately count the vortex 
line elements parallel and antiparallel to the field: 

m m+ m— 

One can assume that the vortices are uniformely distributed in the planes and that the y 
components of their positions are distributed on a line, separated each other by a distance 
AL — -^.Then, the i— th vortex is in the mean position < i?j >= ALi — ^i, with 



_ _N+ N+ 



*~ 2 

By considering N± » 1 one finally finds 

m 

Then (Al) can be written 

Y,t{m) < cos qR{m) >= ^ ' + o{q') 

m. U U 
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VIII. FIGURE CAPTIONS 



Fig.l Some magnetization data in low field, parallel to the c-axis, as a function of temper- 
ature in oriented powders of Yi_ri.CaxBa2Cu30y. The values of the magnetization measured 
from 250 K down to 90 K, with a positive Pauli-like temperature independent term, are not 
reported in the figure. In the insets the blow-up for the estimate of Tc{0) is shown. 

Fig. 2 Constant field (H || c) magnetization vs. temperature in overdoped {x — 0.1 and 
y — 6.96) (a) and in underdoped {x — and y — 6.65) (b) YBCO compounds. For 
comparison in the part a) of the figure the behaviour of Mfi in optimally doped YBCO is 
shown. 

Fig. 3 Isothermal diamagnetic magnetization vs. H, after zero-field cooling (ZFC) to a 
certain temperature above Tc(0). a) sample at x — 0.1, y — 6.96 ( Tc{0) — 73K); b) sample at 
X = 0.2, y = 6.98 ( Te(0) = 49.5K); c) chain disordered underdoped YBCO at y = 6.65 and 
X = { Tc(0) = 62.5K) for ZFC (circle) and field-cooled (FC) (up-triangle) conditions. The 
solid lines in part a) of the Figure correspond to the diamagnetic susceptibility estimated 
in the limit of zero field. The solid lines in Figs. 3a and 3c are the theoretical behaviours 
according the mechanisms described in the text for droplets below and above the local 
irreversibility temperature. 

Fig. 4 Comparison of the magnetization curves Mfi vs. H (after ZFC) in overdoped 
YBCO:Ca {x = 0.1, y = 6.96) with the ones in optimally doped YBCO (inset), for similar 
reduced temperature e. The sohd lines fitting the data in optimally doped YBCO are derived 
from Eq.2 in the text by means of numerical derivation, correspond to anisotropic parameter 
r=0.1 and evidence the 3D linear and 3D non- linear regimes. 

Fig. 5 Susceptibility as a function of the inverse temperature showing the activated tem- 
perature behavior in the sample ai y = 6.96 and x = 0.1. Analogous temperature depen- 
dence has been observed in underdoped YBCO compounds. The dashed lines are obtained 
by transferring above the superconducting temperature T'c(O) the temperature behavior of 
the bulk susceptibihty measured below Tc and by normahzing the data at T ~ 8SK. 
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Fig. 6 Magnetization curves in YBCO:Ca at x — 0.1 and y — 6.96 by cooling in zero 
field to a given temperature [a) T = 75.5K; b) T = 79.5K] above the resistive transi- 
tion temperature and then applying the field (ZFC) and after the application of the field 
at room temperature, cooling at the same temperature and measuring the correspondent 
magnetization (FC). The volume susceptibility in the limit H^O is reported. 

Fig. 7 Relaxation of the raw magnetization after ZFC and then sudden application of a 
magnetic field of 260 G, in YBCO:Ca at x = 0.1 and y = 6.96, at T = 75.5K a) short-term 
relaxation; b) long-term relaxation. From the comparison of the ZFC and FC magnetization 
in H=20G (see inset) an irreversibility temperature of the locally superconducting droplets 
at highest Tc can be estimated around 90 K. In part a) of the Figure the solid line is the 
sketchy behavior of the relaxation of the magnetization detected by Yeshurun et al. (Ref.^^) 
in the critical state, in optimally doped YBCO. 
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IX. TABLE CAPTION 



Table I: Superconducting transition temperature in overdoped and underdoped 
Yi_j;Caa;Ba2Cu30y and estimated number of holes per Cu02 unit 



TABLE I 



X 


y 


UK) 


rih 





6.65 


62.5 


0.12 


0.05 


6.97 


82.0 


0.18 


0.1 


6.96 


73.0 


0.20 


0.1 


6.96 


70.0 


0.21 


0.2 


6.98 


49.5 


0.23 


0.15 


6.10 


34.00 ± 1 


0.07 



i 0.15 6.05 20.00 ±2 0.06 
0.1 6.10 14.00 ±2 0.06 
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